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The behavior of a trapped atomic gas undergoing a Bose-Einstein condensation is strongly influenced by the 
attractive interaction between particles. The role played by such an interaction is analyzed by means of a 
variation method. The critical number of atoms for a metastable state to exist is determined. The aspect ratio 
of condensed bosons in an anisotropic trap, the loss rate of atoms, and the excitation spectrum are studied. The 


vortex state is discussed. [S1050-2947(97)04 104-8] 
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Bose-Einstein condensation (BEC) is a fundamental phe- 
nomenon where a macroscopic number of quantum particles 
is in the ground state of the system at a finite temperature. 
The concept of BEC has been used to explain the phase 
transitions in superfluids and superconductors, and has been 
related to physical phenomena on all scales, including exci- 
tons in condensed matter, Kaons in dense nuclear matter, and 
elementary particles [1]. Although known for almost seventy 
years, BEC was only recently observed directly in dilute 
atomic ®’/Rb, ’Li, and Na vapors [2-4]. These outstanding 
experimental achievements have renewed theoretical interest 
in studying Bose gases with a nonuniform density. 

Among these experiments, perhaps the BEC of ’Li pre- 
sents the greatest challenge to theoretical physics. The 
s-wave scattering length a= —(27.3+0.8)ao, where dg is 
the Bohr radius, of such atoms in the corresponding spin 
state is negative [5], indicating an attractive atom-atom in- 
teraction. It has been claimed that the BEC of attractive 
bosons in free space is impossible [6] because the attraction 
makes the system tend to an ever denser phase. However, for 
a finite number of atoms in a limited volume and a short time 
the situation may be different. For bosons trapped in an ex- 
ternal potential, the potential has an effect that resists the 
tendency to collapse. So, there may exist a metastable BEC 
state under certain conditions. For a harmonic oscillator trap, 
recent numerical results [7—9] for T=0 show that a stable 
condensation can occur for a number of atoms below the 
critical value, which is about 1400 under the condition of 
present experiments. 

Here, we shall study the BEC of a boson gas with an 
attractive interaction at T=0 using a variation method. This 
method was first introduced in Ref. [10] to study the BEC 
ground state in the harmonic trap of a boson system. We 
shall concentrate on the BEC experiment of ’Li, and study 
the stability, aspect ratio, loss rate, and excitation spectrum 
of the condensed state. We also investigate the case when 
vortex lines are present. 

When all atoms are in the same quantum state at T=0, let 
y(r) be the ground state wave function, where 
Sd?r\yr)|?=N. In the mean-field approximation, the 
ground-state energy of the system is given by the Ginzburg- 
Pitaevskii-Gross energy functional [11] 
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where w? and w? are the two angular frequencies associated 
with the external potential of an anisotropic harmonic trap, 
a is the s-wave scattering length, which is negative when the 
atom-atom interaction is attractive, and r =x°+y?, We 
shall consider only the case of a <0. 

The Ginzburg-Pitaevskii-Gross theory is based on the 
pseudopotential form of the atom-atom interaction, and may 
be assumed accurate when the system is dilute. Denote by 
p the boson density. The theory is valid only when 
|a|3p<l1 . This condition is well satisfied for the samples of 
alkali atoms in experiments to date [2—4]. The ground-state 
wave function is determined by minimizing the energy func- 
tional. The variation of the energy functional with respect to 
the wave function gives a nonlinear Schrodinger equation for 
condensed neutral atoms in a harmonic trap. The time- 
dependent [13] and time-independent nonlinear Schrodinger 
equations [12,7] have already been numerically solved. 

The problem now will be solved by the variation with 
respect to a properly chosen trial wave function with adjust- 
able parameters. We first consider the case in the absence of 
vortex lines. Since the atom interaction is rather weak, in the 
BEC state we assume the trial wave function for y in Eq. (1) 
to be [10] 
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where w, and w, are the effective frequencies and are taken 
as variational parameters. Substituting expression (2) into 
Eq. (1), we obtained the ground-state energy 
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(3) 


A physical state corresponds to a stable or metastable point 
of the energy functional. The BEC of a system with an at- 
tractive interaction happens when energy (3) exhibits a mini- 
mal value. When wave functions are restricted to the form of 
the trial function (2), we may write the conditions of a mini- 
mal energy in terms of derivatives of the energy with respect 
to the adjustable variation parameters of the trial function. 
They are 
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By denoting A, =!/a, , A,=o/w, , and A= w/o!) , in- 
troducing the transversal characteristic length of the har- 
monic potential a, =(h/m oy and using expression (3), 
we rewrite the above conditions as 
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Since A, and A, are non-negative, from Eqs. (7) and (8) we 
have 0<A, <1 and 0<A,<1, which means that the trans- 
versal (x-y) and axial (z) root-mean-square radii must be 
smaller than their corresponding characteristic lengths, i.e., 
V(r) <a, and Vz) <(h/mo?) N=q,. 

While Eqs. (7) and (8) are the conditions for stationary 
points, inequality (9) is that for stability. Setting the left-hand 
side of Eq. (9) to zero, i.e., replacing the inequality by an 
equation, we can find the critical conditions for BEC. Taking 
into consideration the experimental parameters w? =163 Hz 
and o= 117 Hz of the harmonic trap in Ref. [3], we obtain 
X~0.71779 and N it~ 0.2808 V2 ma, /|a|~ 1444. The BEC 
happens only when the number N of atoms is smaller than 
this critical value. 

In experiments the onset of BEC is signaled by a narrow 
peak centered at zero velocity on the top of a broad thermal 
velocity distribution [2—4]. For the case of positive scattering 
length the peak exhibits a nonthermal, anisotropic velocity 
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FIG. 1. Root-mean-square radii along the x axis (dashed line), 
along the z axis (dotted line), and aspect ratio R (solid line) for 
different numbers of atoms. The lengths are in units of a, . 


distribution. The parameter to characterize the asymmetry of 
the velocity distribution function is the so-called aspect ratio 
R= Vp?) (p2). The momentum distribution of particles in 
the trap is given by f(p)=|fd?rexp(—ip-r/h) W(r)|7. Our 
Gaussian trial wave function (2) then leads to the Gaussian 
momentum distribution f (pP) ~expl- p/w u + lo,)], with 
(p2)=(pi)=mhw,/2 and (p?)=mhw-/2. Finally, we find 
the aspect ratio R= w,/@, = VMMA, /A,. The root-mean- 
square radii VaI a], VII a, and aspect radio R as func- 
tions of the number of atoms are shown in Fig. 1. We see 
that when the number of atoms increases, the attractive in- 
teraction between atoms becomes important, the root-mean- 
square radii decrease, and the aspect radio approaches 1. 

It has been pointed out that the transition to an unstable 
state, where the kinetic energy of ’Li atoms in a harmonic 
potential can no longer stabilize the wave function and pre- 
vent the atoms from collapsing, may occur due to a quantum 
tunneling [14]. This tunneling is analogous to the ordinary 
quantum tunneling of a particle in an unstable potential, but 
is much slower [15]. Neglecting the tunneling and consider- 
ing only the loss rates due to the two-body dipolar collision 
and three-body recombination collision, we have the total 
loss rate [16,17] 


: (10) 
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where K is the two-body dipolar loss rate coefficient, and 
L the three-body recombination loss rate coefficient. Substi- 
tuting the wave function (2) into Eq. (10), we find 
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Here we have neglected the interaction of the condensed at- 
oms with the background thermal gas. This y gives a lower 
bound to the loss rate. The loss rate for different numbers of 
atoms is shown in Fig. 2. The parameters used for the figure 
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FIG. 2. Loss rate of condensed atoms as a function of the num- 
ber of atoms. 


are K=1.2X 107! cmêîs™! and L=2.6X 10778 cmfs™! as 
in Ref. [9]. We see that the loss rate increases very rapidly 
with an increasing number of atoms, especially near N erit- 
In order to obtain some analytic results, we consider the 
isotropic harmonic oscillator potential by setting 


w? = w= wp. Equation (3) then reduces to 
L z 0 q 
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Correspondingly, the conditions for the energy function to be 
minimal are 
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which determine the critical number of trapped atoms to be 
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where rọo= Val mo. In order to compare with the case of an 
anisotropic potential, we take wọ = (20? + w?)/3. For the ex- 
perimental parameters of Ref. [3], N or™ 1446. 

From Eq. (13), we have 
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Taking inequality (14) into account, we see that if the inter- 
action is attractive, i.e., a <0, then dw/dN>0. That is, the 
root-mean-square radius decreases with an increasing num- 
ber of atoms. 

Based on the ground state found above we now discuss 
the elementary excitation spectrum. According to the Bogo- 
liubov approximation [18], the effective energy functional of 
an excited state is 
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where u(r) and v(r) describe the components of the elemen- 
tary excitation wave function. The commutation relation re- 
quires 


| dr[u*(r)u(r)—v*(r)u(r)]=1. (18) 


For simplicity, we study only the case of the spheric har- 
monic trap. We choose the trial wave functions of the exci- 
tation components in the form of the spheric harmonic oscil- 
lator wave function Vn, dim With quantum numbers 
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where œ is an adjustable scaling factor of variation. For 


(n,,l,m)=(0,1,0), we have 
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The energy functional (17) reduces to 
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where Apio= W9/Wo10; A; = wo/w, and w010 7 whim. Simi- 
larly, for (n,,/,m)=(1,0 0), we have 
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where 
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FIG. 3. Excitation frequencies in units of wọ of the modes 
(0,1,0) and (1,0,0) for different numbers of atoms. 
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A, is given by Eq. (13). Minimizing the energies of Eqs. (21) 
and (23) with respect to the variation parameters v, Apgjo, 
and Ajgg, we can determine the excitation frequencies for 
the modes (0,1,0) and (1,0,0). The results are shown in Fig. 
3, from which we see that while for the mode (0,1,0) the 
excitation frequency increases slightly with an increasing 
number of atoms, for the mode (1,0,0) the excitation fre- 
quency decreases. This is in agreement with Ref. [9]. 

We now turn to a discussion of vortex states. By assum- 
ing that all the atoms rotate around the z axis with a quan- 
tized circulation, the wave function of the condensed vortex 
state can be written as [19] 


W(r) = p(r)exp[i® (1) ], (25) 


where o(r) = Vem) is the modulus, and the phase ®(r) acts 
as a velocity potential, i.e., V=(/m)V®(r). We further 
take B(r)=«d, where ¢ is the angle around the z axis, and 
k the integer quantum number characterizing the circulation. 
The angular momentum L, along the z axis is N«h. For the 
wave function of this form, the Ginzburg-Pitaevskii-Gross 
energy functional (1) becomes 
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For simplicity, we shall consider only the vortex state with 
k=1. We choose the trial wave function to be of the same 
form as the wave function for noninteracting atoms given in 
Ref. [7] 
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FIG. 4. Critical angular velocity in units of w? for the vortex 
state of k=1 in Li vapor as a function of the number N of atoms. 
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For a higher « a different trial wave function g(r) should be 
taken. Substituting expression (27) into (26) we obtain 
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By similarity to the case in the absence of vortex lines, the 
conditions for the energy function to be minimal are 
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From these conditions we find that for the given experimen- 
tal parameters, when N|a|/./2 77a, >0.8457 the energy func- 
tion cannot be minimal. That is, when the number of atoms 
N>N., vo™ 4349 there is no stable or metastable vortex state 
with k=1. This agrees with Ref. [7]. Compared with the 
condensed state without vortices, the critical number of 
trapped atoms in the vortex state increases. When all the 
atoms rotate around the z axis, the Coriolis force is favorable 
for preventing attractive atoms from collapsing. 

Furthermore, the critical angular velocity Q of rotation 
can be obtained from 6E—L,0,<O [20], where ôE is the 
energy gap between the ground and vortex state, i.e., 
E=E —E ground: The critical value is 


vortex 
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For a given number of atoms, A, and A, can be calculated, 
and from Eqs. (3) and (28) the ground-state energy and 
vortex-state energy can be further computed. The critical an- 
gular velocity for the vortex state of k=1 in ’Li vapor as a 


function of the number of atoms is shown in Fig. 4, from 
which we see that the critical angular velocity, or vortex 
frequency, increases as the number of atoms grows. This is 
different from the case when the interaction between atoms 
is repulsive [7,21]. 
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